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Abstract 

We extract diffractive parton densities from diffractive, deep inelastic (DIS) ep data from the 
ZEUS experiment. Then we use these fits to predict the diffractive production of jets and of 
JU’s and ZA in pp collisions at the Tevatron. Although the DIS data require a hard quark 
density in the pomeron, we find fairly low rates for the Tevatron processes (a few percent of 
the inclusive cross section). This results from the combined effects of Q 2 evolution and of 
a normalization of the parton densities to the data. The calculated rates for W production 
are generally consistent with the preliminary data from the Tevatron. However, the jet 
data from CDF with a “Roman pot” trigger are substantially lower than the results of our 
calculations; if confirmed, this would signal a breakdown of hard-scattering factorization. 
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1 Introduction 


In view of counterexamples |l| to the conjecture of factorization || of hard processes in 
diffractive scattering, it is important to test [|3J factorization experimentally. In this paper, 
we present some results to this end. Specifically, we present some preliminary fits to data 
on diffractive deep inelastic scattering [[| , and use these fits to predict cross sections in hard 
diffractive processes in pp collisions, with the assumption of factorization. 

We recall that diffractive events are characterized by a large rapidity gap, a region in 
rapidity where no particles are produced. We are concerned with the case where there is 
a hard scattering and the gap occurs between the hard scattering and one of the beam 
remnants. Such hard diffractive events have been observed in deep inelastic scattering (DIS) 
experiments ||, and are found to have a large rate: around 10% of the inclusive cross section. 
Diffractive jet production in pp collisions was earlier reported by the UA8 collaboration [ 6 j, 
but under somewhat different kinematic conditions (larger |f|[j). There was also a report 
of diffractive bottom production [0. Now, more diffractive data are being gathered from a 
variety of lepto-hadronic j|, §, j 9 j and hadronic processes m 0 0 0 0 but with 
substantially smaller fractions in the case of the diffractive production of jets and weak vector 
bosons in pp interactions. 

The QCD-based model that we use in our calculations is the one due to Ingelman and 
Schlcin []2J, where diffractive scattering is attributed to the exchange of a pomeron — a 
colorless object with vacuum quantum numbers. The pomeron is treated like a real particle, 
and so a diffractive electron-proton collision is considered to be due to an electron-pomeron 
collision. Thus hard cross sections involve a hard scattering coefficient (or Wilson coefficient), 
a known pomeron-proton coupling, and parton densities in the pomeron. Similar remarks 
apply to diffractive hadron-hadron scattering. 

The parton densities in the pomeron can be extracted from diffractive DIS (F 2 ) mea¬ 
surements. Since the pomeron is isosinglet and is its own charge conjugate, there is only a 
single light quark density to measure; one does not have the complications of separating the 
different flavors of quark that one has in the case of the proton. Scaling violations enable 
one to determine the gluon density. The HI collaboration has already presented |J a fit of 
this kind. This type of data sufficiently determines the quark but, at present, it only weakly 
constrains the gluon density in the pomeron. The first experimental evidence for the gluon 
content of the pomeron was found by the ZEUS collaboration by combining their results on 
the diffractive structure function in deep inelastic scattering JfJ and their measurements of 
diffractive inclusive jet photoproduction || under the assumption of factorization between 
these two processes. In this paper we use each of several different fits that we have made to 
diffractive F 2 data 0 from ZEUS. 

Our fits to the DIS data are made with full NLO calculations, but, given the level 
of accuracy that is needed at present, we only use lowest order QCD calculations for the 
hadronic processes. In the past, Ingelman and Schlein || and Bruni and Ingelman [[HJ have 
made similar calculations for one of the hadron-induced processes that we consider here. 
Their results have provided a commonly used benchmark in the phenomenology of these 
processes. They provide a choice of either ‘hard’ or ‘soft’ distributions of partons in the 


1 By t we mean the invariant momentum-transfer-squared from the diffracted hadron. 
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pomeron, according to the (3 —> 1 behavior.^ The hard distributions give larger diffractive 
cross sections. At that time, there were no data to determine the distributions. We will 
find that although the quark distributions preferred by the DIS data are hard, our cross 
sections are substantially below those predicted by Bruni and Ingelman. As we will see, the 
lower cross sections occur for several reasons, particularly: a correct normalization of the 
distributions to DIS data and the incorporation of Q 2 3 evolution. 

How well our predictions match up with the data in hadron-hadron collisions will be a 
statement on the validity of factorization of the diffractive hadronic cross sections. If there 
is good agreement, comparison of our results with measured cross sections will also provide 
a good test of the pomeron parton distribution fits from diffractive DIS. In addition, the 
results from diffractive jet photoproduction |J will provide greatly improved constraints 
on the gluon content of the pomeron, and will substantially improve the accuracy of our 
predictions.^ 

This paper is organized as follows. In section we present the formulae used to calculate 
the various cross sections. We also discuss the kinematics and phase space cuts that were 
used. In section |3|. we present our fits to the diffractive deep inelastic data. Then in sections 
[| and ^ we present and discuss the results obtained for vector boson production and jet 
production, respectively. We will find that production rates for these diffractive events are 
no more than a few percent, at most, of the nondiffractive ones. Finally, we summarize our 
findings in section |j. 

Other fits to the diffractive structure functions measured by HI have been made by 
Gehrmann and Stirling |l8| and by Kunszt and Stirling |T9|. Golec-Biernat and Kwiecinski 
PD|| assumed a parameterization of the parton densities in the pomeron and found it to be 
compatible with the HI data on diffractive DIS. Their quark densities are about 30% smaller 
than ours, and they required the momentum sum rule to be valid. The new features of our 
work are a fit to the ZEUS data, and a calculation of the cross sections for diffractive jet 
and W and Z production. 


2 Kinematics and Cross Sections 

The diffractive processes that we consider here are the production of W and Z bosons and 
of jets in pp collisions. In addition we consider W production with explicit consideration of 
the distribution of the hnal state leptons. Schematically, these are 

p(Pi) +p{p- 2 ) —> (W or Z)+p + X, 
p(Pi)+p(p 2 ) —> jet +p + X 

p(Pi)+p(p 2 ) —> (W —> l + v) +p + X. (1) 

2 Here, (3 is the fraction of the pomeron’s momentum that is carried by the struck parton. 

3 As explained in Refs. §!• we consider it more likely that factorization is valid in diffractive deep 
inelastic and direct photoproduction than in diffractive hadron-induced processes. Therefore, we prefer 
to determine the parton densities in the pomeron from deep inelastic scattering and from direct photo¬ 
production, and then to treat the hadron-induced diffractive processes as providing tests of factorization. 
Goulianos’ proposal 0 to renormalize the pomeron flux in an energy dependent way could be regarded 
already as evidence that factorization is likely to break down. 
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We take the pomeron to be emitted from the antiproton and the positive z-axis to be along 
the antiproton’s direction. The center of mass energy is y/s, where s = (p\ + P 2 ) 2 - 


2.1 Diffractive jet production 

Consider the diffractive cross section for the production of a jet with rapidity y, in a hadron- 
hadron collision. We will assume hard-scattering factorization i a- The lowest order 
hard-scattering process is 2 —> 2 at the parton level, and results in a cross section of the 
form 

/ 2/max 

dy 

-'min 

( 2 ) 

where the sum is over all the active parton (quark, antiquark and gluon) flavors. The 
integration variables are Et, the transverse energy of the jet, y\ the rapidity of the other jet, 
and xp, the momentum fraction of the pomeron. The momentum fractions of the partons, 
relative to their parent proton and pomeron are 



dxp fp/pi^Xp, p) fa/pfaa, h) fb/w(x b , !l)X a X b 


ddi 


et 

ab 


dt 


X a 

X b 

The functions f a / P (x a ) and f b / v (x b ) 
pomeron, respectively, while fp/p is the “flux of pomerons in the (anti)proton”, to be dis¬ 
cussed below, da^/dt is the partonic hard scattering coefficient and /i is the factorization 
scale, which we set equal to E r . The specific limits used for the integral in xp, as well as 
those for the rapidity y' and transverse energy Et will be given later. 

The diffractive cross section given by Eq. (|2|) has the same structure as the factorized 
form of the corresponding nondiffractive cross section, except for the pomeron flux factor 
and the parton densities in the pomeron. The same hard scattering coefficient and nucleon 
parton distribution functions appear in both cross sections. 

The pomeron flux factor, fp/p, is related to the pomeron-proton coupling measured in 
proton-proton elastic scattering. The t dependence (which we integrate over) and the xp 
dependence are thereby determined, but there is a controversy as to the (constant) nor¬ 
malization needed to treat the exchanged pomeron as if it were a particle. Since the same 
normalization factor appears in all our cross sections, its value will be irrelevant to our phe¬ 
nomenology. Any change in the normalization factor is completely compensated by changing 
the parton densities by an inverse factor, and we obtain the parton densities from fitting a 
set of data without any a priori expectations as to their normalization. 

However, the normalization does affect the question of whether the momentum sum rule 
is obeyed by the parton densities in the pomeron. Since it is not at present understood 
whether the sum rule is a theorem, this issue will not affect us. The momentum sum rule is 
not assumed in any of our fits. 

4 These are number densities, not momentum densities. 
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{e y + e y '). 


E T 
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are the distributions^] of partons in the proton and 
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There are two pomeron flux factors that are commonly used, Ingelman-Schlein (IS) []2| 
and Donnachie-Landshoff (DL) Since the parton densities will be obtained using the 
DL factor, we use the same factor in predicting other cross sections. The DL flux factor is 
given by 


/p+ Op) = 


r dt ?%\ 

■4m 2 — 2.8 1 

( 1 \ 2_ 

/-oo 4 P 1 

. 4 m^ — t 

a-t/o.7/ J 


x 


1-2 a(t) 


( 4 ) 


where m p is the proton mass, /3q ~ 1.8 GeV -1 is the pomeron-quark coupling and a(t) = 
1.085 + 0.25t is the pomeron trajectory. The integral is over the invariant momentum transfer 
carried by the pomeron, since in the present generation of measurements, only the value of 
the longitudinal momentum of the pomeron is measured. Since the t distribution is steeply 
falling, only values of \t\ under about a GeV 2 are significant. 

The cross section given by Eq. (|2;) has contributions from a range of subprocesses. The 
indices a, b labeling the incoming partons range over the gluon and all the flavors of quarks 
and antiquarks. The lowest order form of the partonic cross section da^ b /di may be found 
in |22| . 


2.2 Diffractive W and Z production 

The cross section for the diffractive production of weak vector bosons, is given by 


cr 


w,z 



dx b f ¥/p (x v )f b/v (x b )f a/p {x a )C 2 ab 


sx b x pL 






VB 


( 5 ) 


where x b is now the momentum fraction of the parton from the pomeron, so that the mo¬ 
mentum fraction of the other parton (in the proton) is x a = My B /x b XpS. The minimum 
value of x b is = My B /xps. Also, M VB = M w or M z is the vector boson mass, and Gp 
is the Fermi constant. For W bosons, C')|' = V ab , the relevant Cabibbo-Kobayashi-Maskawa 
matrix element, while for the Z boson, 



2\e b \ sin 2 0 W + 41easin' 1 6 W , 


( 6 ) 


where e b is the fractional charge of parton b and 0\y is the Weinberg or weak-mixing angle. 


2.3 Diffractive production of leptons from the W 

Since leptonic decays of W bosons include an unobserved neutrino, it is useful to compute 
the distribution of the observed charged lepton. The general formula for the distribution of 
leptons from W production has the same form as that for jet production, Eq. (|J). In the 
case of W production, at lowest order there are no gluon contributions. 

For the specific process p + p —► ( W~ —> e + u e ) + p + X, we have 

,/D e P 2 

* -dTr^V 2 b 5{x a x b s - M 2 w )u 2 , (7) 

at OM\yL w 
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where Mwi^w) is the mass (width) of the W boson, V ab is the Cabibbo-Kobayashi-Maskawa 
matrix element and 

u = — XbXpy/sE T e~ y . ( 8 ) 

We used the narrow width approximation in Eq. (0). Using Eq. (|7|) in Eq. (^|), one obtains 


da lep 

dy 




dErfr/pM fb/p(x b )f a/p (x a )V 2 b 


u 2 G 2 F 

6sT\y[(M\y / 2Et) 2 


where x a , x b are now given by 



1 


(9) 


X a 


Xb 


M w e y M w j / M w \ 2 
y/s 2 E t V \2Et) 

M 2 W 1 
s x a x P ' 


( 10 ) 


We have suppressed the scale dependence of the functions fa/j in Eqs. (f|) and (|9|); we set 
the scale equal to the vector boson mass. A similar equation may be obtained for the W + 
cross section. 


2.4 Inclusive cross sections 


Since we are particularly interested in the percentage of events that are diffractive, we also 
need to calculate the inclusive cross sections, that is, the ones without the diffractive re¬ 
quirement on the final state. The analog to Eq. (^) for the inclusive cross section for jet 
production is the standard formula 


da J et > 


incl 


dy 


E 


dE T 2E 7 


dy fa/p^Xai y) fb/p(x b i y)x a x b 


dot 


jet 

ab 


dt 


( 11 ) 


where x a is given in Eq. @ while x b is now x b = ( e y + e y ')Er j\fs. 
For the leptons from W~ production, the inclusive version of Eq. 


is 


d(J lep, incl 

dy 


E 


/ zpmin 
a, b J 


dEx fb/p{x b: fl) fa/p(x a , fi)V ab 


u 2 G 2 F 


6sT w [(M w / 2E t ) 2 - l] 1 / 2 


( 12 ) 


with a similar equation for W + production. In Eq. (|12D , u = — x b yfsE T e~ y ) x a is as defined 
in Eq. ([H]) while x b is now given by x b = Myy/x a s. 

The analog to Eq. @ for the inclusive total cross section for vector boson production is 


<T W ’ Z = 


'y ^ I dx b f b / p (x b )f a / p (x a )C~ b 


a.b 


SX 6 L 


V2| G , 


■M, 


VB 


( 13 ) 


where C ab and Mvb are as defined above, = My B /s and x a = My B /x b s. 
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3 Partons in the Pomeron 


We have made five fits of parton densities in the pomeron from the ZEUS data in Ref. Q 
for the diffractive structure function F 2 . Diffractive structure functions are related to the 
differential cross section for the process e + p—*e + p + X\ 


d C^diff 

df3dQ 2 dxpdt 


2n a 2 

~w 


1 + (1 -y) 2 ]F? {A) -y 2 F?V>}, 


(14) 


where corrections due to Z° exchange and due to radiative corrections have been ignored. 
Here xp and t are the same as in the previous sections, Q 2 and y are the usual DIS variables, 
and (3 = x^/xp, with x bj being the usual Bjorken scaling variable of DIS. The diffractive 
structure function F 2 < ' 4 \/3,Q 2 ,xp,t) is assumed to obey Regge factorization, so that it is 
written as a pomeron flux factor times a pomeron structure function: 

F 2 Dii \p,Q 2 ,x T ,t) = Mx r ,t)F'(P,Q 2 ). (15) 


The momentum transfer t is not measured, so the data are actually for the structure function 
integrated over t: 

F 2 d{3 \(3,Q 2 ,xp) = [° dtF^\/3,Q 2 ,xp,t). (16) 

Furthermore, the actual fits are to F 2 ((3, Q 2 ), which is obtained by integrating F 2 ^ ((3, Q 2 , xp) 
over the measured range of xp, 6.3 • 1CT 4 < xp < 1CT 2 , using the fitted xp dependence 0. 
That procedure, as noted in [[J, assumes that a universal xp dependence holds in all regions 
of (3 and Q 2 . 

Hard scattering factorization gives F 2 in terms of parton densities and hard scattering 
coefficients in the usual fashion: 


F 2 ((3, Q 2 ) = e 2 a pf a/ p(p) + NLO corrections. (17) 

a 

Since the outgoing proton is not detected, the data include contributions where the proton 
is excited to a state that escapes down the beam-pipe and thus misses the detector. Excited 
states up to about 4 GeV pass the diffraction selection cuts, and the experimenters estimate 
that there is a contribution of (15 ± 10)% to the measured diffractive F 2 from such “double- 
dissociative” events. This point is significant when we compare predictions obtained using 
our fits to data where the diffracted proton is detected, as in Sect. [|. 

Each of our fits is represented by a parameterization of the initial distributions at Qq = 
4 GeV 2 for the u, u, d, and d quarks and for the gluon. The other quark distributions 
are assumed to be zero at this scale. The fits were made with NLO calculations (with full 
evolution and with the number of flavors set equal to 4) and with the pomeron flux factor 
chosen to be that of Donnachie and Landshoff|] The program used to perform the evolution 
was that of CTEQ [p3|1 - 

5 The flux factor is a common factor in all the cross sections we compute, and we do not assume a 
momentum sum rule for the parton densities in the pomeron. Therefore, the choice of flux factor does not 
affect our predictions for the hadron-induced processes provided only that the pomeron trajectory function 
a(t) is correct. 
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NLO —QCD fits to Fj ZEUS Data 



O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


(3 

Figure 1: The (3 dependence of the diffractive structure function F 2 D measured by ZEUS, together 
with our fits. The inner, thick error bars represent the statistical errors on the data, and the 
outer, thin error bars represent the systematic and statistical errors added in quadrature. Fit A 
is represented by the dashed line, fit B by the thin solid line, fit C by the dot-dashed line, fit 
D by the dotted line, and fit SG by the thick solid line. Note that fits B and D are essentially 
indistinguishable. 


Four of the fits, labeled “A”, “B”, “C” and “D”, use conventional shapes for the initial 
distributions. The final fit has a gluon distribution that is peaked near (3 — 1, as suggested 
by the fit [§] exhibited by the HI collaboration; we call this our “singular gluon” fit, SG. We 
show our fits in Figs. [T] and [2]. 

Our first fit, A, has as its initial distribution a hard quark distribution (proportional to 
[3(1 — /3)) with no glue, and hence one adjustable parameter, whose value is the result of the 
fit: 


= 0.585/3(1 -/3), 

PftlMQl) = 0. (18) 

This fit gives the dashed line in the figures. It represents the shape of the diffractive structure 
function F 2 D moderately well, but there are noticeable deviations, both as a function of Q 2 
and of (3. 

For the Q 2 -dependence, Fig. ||. we see that fit A results in an F 2 D that at large (3 decreases 
with Q 2 , whereas the data show a tendency to increase. At small /3, both the data and the 
fit rise with Q 2 , but the data rise more rapidly. This suggests adding in an initial gluon 
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NLO —QCD fits to F° ZEUS Data 



1 O 20 30 40 50 60 


Q 2 [GeV 2 ] 

Figure 2: The Q 2 dependence of the diffractive structure function F 2 D measured by ZEUS, 
together with our fits. The description of the lines and error bars is the same as in Fig. jT| 
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distribution, whose effect is to make the quark distributions rise. As to the /3-dependence, 
we see from Fig. jl], that the fit gives an F 2 D that is rather higher than the data at large /3 
and rather lower at small /3. This suggests that an admixture of soft quarks, more peaked 
at small (3 would improve the fit noticeably. 

So we try each of these additions in turn, and then together. First, in fit B we allow a 
hard gluon in the initial distribution, with the result: 

Pfq/A^Qo) = 0.516/3(1-/3), 

= 12.28/3(1-/3). (19) 

Since the gluons only affect deep inelastic scattering in next-to-leading order (NLO), a very 
large gluon distribution is needed to produce a substantial effect on F 2 D . The resulting Q 2 
evolution has a much more satisfactory shape, as can be seen from the thin solid line in 
Fig. |2j. In addition, the strong gluon distribution has the effect of biasing the evolved quark 
distribution towards smaller /3. Thus the shape of idP as a function of /3 is also improved, 
Fig. HQ 

Next, in fit C, represented by the dot-dashed line, we see the effect of adding a soft term 
to the quark distribution without allowing an initial gluon distribution. In general we expect 
a soft quark term, since Regge theory predicts that a quark distribution at small /3 behaves 
approximately as /3/ ~ 1. We find the best fit of this kind to be: 

U {')-.(= 0.470/3(1 -/3) + 0.080 

Pf°lr(P,Ql) = 0 . ( 20 ) 

The normalization of the hard quark term is reduced by about 20%, compared with fit A, 
and the added soft term clearly improves the shape of the /3 distributions, in Fig. |l]. It should 
be clear that at moderate and large values of /3, above about 0.2, there is a dominant hard 
quark term. Moreover a soft term, at small /3 is needed, no matter whether it is intrinsic to 
the quark distribution or whether it is generated dynamically, from evolution controlled by 
the gluon distribution. As is to be expected, the soft quark term does nothing to improve 
the Q 2 dependence in Fig. §. 

We next try both a soft quark term and a hard gluon term. This results in fit D (shown 
as the dotted line in the figures): 

= 0.512/3(1 -/3) + 0.005 (1 -/3) 2 , 

Pf?,AP,Ql) = 11-65/3(1-/3). (21) 

This is indistinguishable from fit B, particularly given the errors. The important fact is that 
if the pomeron were quark dominated, QCD predicts that F 2 D decreases with Q 2 at large /3. 
The only way of undoing, or even reversing, this evolution is to have a large initial amount 
of glue. No extra soft quark term is required by the data. 

Given the number of points and their errors, we can conclude that 

• Hard quark and gluon distributions are preferred. 

6 Remember that the initial Q 2 is Q g = 4GeV 2 , so that evolution already has a significant effect at the 
lowest value of Q 2 for which we use data. 
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• The normalization of the quark distribution is reasonably well determined. 

• Probably a substantial gluon distribution is preferred. 

However, even the pure quark fit, A, gives a satisfactory y 2 — see Eq. ( j23| ) below. It should 
be remembered, however, that y 2 is not the only measure of goodness-of-fit. On the other 
hand, the error bars in the lowest plot in Fig. |2] indicate that reduction of the systematic 
errors in the measurement of OF®/dQ 2 is urgently needed to improve the determination of 
the gluon distribution. 

Finally, we recall that the HI collaboration has shown fits to their data with a gluon 
distribution that is very peaked close to (3 = 1. So we have tried a similar fit. The result, fit 
SG (shown as the thick solid line in the figures), is: 

(3fq/v(Pi Qo) = 0.354/3(1 — (3), 

Pfg/ G AP,Ql) = 70.756/3 8 (1 — /3) 0 ' 3 . (22) 

The exponents in the gluon density of Eq. (|22|) were chosen to try and match the singular 
gluon parameterization derived by HI the ZEUS data, at least at this time, are not able 
to determine definite values for these exponents. To avoid possible numerical problems with 
the evolution code we have chosen a shape for the gluon distribution that has a peak close 
to (3 — 1 but that is not actually singular there. 

The singular fit results in an improved y 2 — see Eq. (|3]) below — presumably because 
the resulting F 2 D increases more rapidly with Q 2 at large /?. Given the size of the systematic 
error bars in the relevant plot, the bottom one of Fig. |2], this reinforces the need for improving 
the systematic errors on dF® / dQ 2 . 

In each of our fits, the parameters of the fit are the nonzero coefficients, with the expo¬ 
nents being held fixed. The parameters are determined by minimizing the y 2 between the 
theory and the data. The systematic and statistical errors have been added in quadrature, 
although that is not an ideal procedure. All the fits give a good y 2 and the fits with the 
large gluon are somewhat preferred, particularly the singular gluon fit: 



Fit A 

Fit B 

Fit C 

Fit D 

Fit SG 

y 2 /d.o.f 
Statistical y 2 /d.o.f 

5.7/10 

18/10 

1.7/9 

3.8/9 

2.5/9 

6.0/9 

1.7/8 

3.8/8 

1.2/9 

3.1/9 


The fact that the y 2 per degree of freedom for all these fits is much less than unity presumably 
indicates that the systematic errors have not been treated correctly, for example, as regards 
point-to-point correlations. However, even if only the statistical errors are taken into account, 
4 of the 5 fits have a worryingly low y 2 per degree of freedom, as can be seen from the last 
line. There is only a few per cent probability of getting fits this good, even if the theory is 
exactly correct. We take this as an indication that one should look more carefully at the 
computation of the errors on the structure functions. 

Independently of the quantification of the errors, the most important feature from our 
point of view is that the fits give rather stable normalizations to the quark distributions. The 
second important feature is that a very large gluon distribution is preferred, perhaps strongly 


11 








peaked at large (5. This is in agreement with the results presented by the HI collaboration 
||. The momentum sums Z a fo dflPfa/ p(/3) are given in the following table: 



Fit A 

Fit B 

Fit C 

Fit D 

Fit SG 

Quarks 

0.39 

0.34 

0.42 

0.35 

0.24 

Gluons 

0 

2.05 

0 

1.94 

3.57 

Total 

0.39 

2.39 

0.42 

2.29 

3.81 


Although the allowed gluon distributions cover a wide range, the overall normalization of 
the quarks is determined to about ±10%, except in the case of the singular gluon fit, where 
the quarks are brought down by about 30%. 


4 Numerical Calculations of W and Z Production 


For the calculations in this section, the factorization scale in the parton distributions was 
set to M vb . The values of the electroweak parameters which appear in the various formulae 
were taken from the particle data handbook [EM|, and we use only four flavors (u, d, s, c) in 
the weak mixing matrix, with the Cabibbo angle 6% = 0.2269. 


4.1 Comparison to previous calculations 


Bruni and Ingelman [|T6[ computed diffractive W/Z cross sections up to 0(a s ), i.e., including 
gluon contributions. These calculations neglected any Q 2 evolution of the parton distribu¬ 
tions in the pomeron. At i/s = 1800 GeV, they obtained the following diffractive fractions 
(R = cr dlff /cr mcl ): R w + +w - ~ 17% and R z — 15% for total W and Z production, respec¬ 
tively. These rates are substantially larger than the few percent quoted in preliminary CDF 
results H 
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As we will now explain, when one uses evolved pomeron parton densities from the above 
fits to data from the ZEUS experiment, one obtains substantially smaller rates than the 
Bruni-Ingelman ones. To understand these small rates, we first verify that we can reproduce 
the Bruni-Ingelman results. For these we used their unevolved hard quark distribution in 
the pomeron (given by their Eq. (4)), the same cut on x^-. = 0.1, the EHLQ1 parton 

distributions in the proton and the Ingelman-Schlein (IS) flux factor:[] 

1 /_ - « - ( 25 ) 


fv/ P fo)= /^^(6.38e 8t ±0.424e 3t ). 


Next, we evolved their pomeron parton distributions and recalculated the cross sections. Fi¬ 
nally, to provide our best estimates of the rates, we repeated the calculations using CTEQ3M 
for the parton densities in the proton/antiproton and using our fits for the parton densities 
in the pomeron, all with proper evolution.^] All the inclusive cross sections were calculated 
using Eqs. ([5]) and (|T5| ) . The results we obtained are summarized in Tables |T] to |3j. 

' Note that since our purpose in using the IS flux is to compare our results with the Bruni-Ingelman 
calculations, we have used a pomeron intercept of unity instead of the more accurate value used in the DL 
factor Eq. (|J). 

8 We evolved the BI distribution from Q\ = 5 GeV 2 (as with the EHLQ1 distributions), while our fits 
were evolved from Qq = 4 GeV 2 . 
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Ref. |16) LO 

EHLQ1 

EHLQ1 

CTEQ3M 

w + + w~ 

14000 

14332 

18150 

z 

4400 

4407 

5383 


Table 1: Inclusive cross sections a w,Zmd (pb) for weak vector boson production. 


Pomeron: 

Proton: 

BI IB 

unevolved 

EHLQ1 

Our BI 
unevolved 
EHLQ1 

BI 

evolved 

EHLQ1 

Fit A 
evolved 
CTEQ3M 

Fit D 
evolved 
CTEQ3M 

W+ + W~ 

2800 

2768 

2025 

518 

844 

z 

760 

738 

520 

133 

204 


Table 2: Diffractive cross section <j w ’ ZdlS (pb) for weak vector boson production, with x™ ax = 
0 . 1 . 


First, in Table [3], we show the inclusiv cross section, <r mc , which will give the denomina¬ 
tor for the fraction of the cross section which is diffractive. Our results are shown in the last 
column. These are 20% to 30% higher than those of Bruni and Ingelman (column 2). We 
have verified (column 3) that this increase is due to the use of the more up-to-date CTEQ3M 
densities in the proton instead of the EHLQ1 densities used by Bruni and Ingelman. Given 
the current accuracy of the diffractive data, we did not bother to make NLO calculations; 
generally the effect would be to increase the cross sections by some tens of percent. 

The diffractive cross sections cr w,ZdlS are shown in Tables and ||. In the columns 
labeled ‘BE, we used the Bruni-Ingelman parton density in the pomeron and the EHLQ1 
parton densities in the proton. In the other columns we used our fits for the parton densities 
in the pomeron together with the CTEQ3M parton distributions in the proton. First, we 
use the same cut x™ ax = 0.1 that was used by Bruni and Ingelman, to produce Table |2j. 
However, this allows Xp to be rather larger than where pomeron exchange is expected to 
dominate. So we also made calculations with ;rp ax = 0.01, for which the results are shown 
in Table |3]. 

In column 3 of Table we show our results when we use the same unevolved parton 
densities as Bruni and Ingelman; we agree with their cross sections. Then we repeat the 
calculations but with correctly evolved parton densities in the pomeron, with the Bruni- 
Ingelman formula being used as the initial data for the evolution at Qq = 5GeV 2 . We 
see that this leads to about a 30% reduction in the cross section. The diffractive fraction 
obtained from the evolved BI pomeron parton distribution, using column 3 of Table |T] for 
fjW,z mci^ j g abQuf; 14% for W production, compared with the 19% that is obtained using the 
unevolved BI pomeron distributions. The corresponding percentages for Z production are a 
little smaller: 12% and 17%. 

In the last two columns of Tables |2] and [3] we present the results when two of our fits to 
ZEUS data are used. Fit A is the one with a simple hard quark distribution and no glue as 

9 I.e., diffractive plus non-diffractive. 
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Pomeron: 

Proton: 

BI (evolved) 
EHLQ1 

Fit A (evolved) 
CTEQ3M 

Fit D (evolved) 
CTEQ3M 

W+ + W~ 

52.3 

12.8 

13.9 

z 

6.6 

1.6 

1.6 


Table 3: Diffractive cross section cr w,ZdlS (pb) for weak vector boson production, but now with 
xjT* = 0.01. 



BI (unevolved) 

= 0.1 

Fit A 

xj? ax = 0.1 

Fit D 

= 0.1 

Fit A 
= 0.01 

Fit D 

X ” ax = 0.01 

W+ + W~ 

19% 

2.9% 

4.7% 

0.07% 

0.08% 

z 

17% 

2.5% 

3.8% 

0.03% 

0.03% 


Table 4: Diffractive fractions 


the initial values, while fit D, together with fit B, has the smallest y 2 ■ Although this latter 
fit has a very large gluon distribution which carries about seven times as much momentum 
as the quarks, the preference for a large amount of glue is mild—see the last two columns of 
([231) . Furthermore, at Qq = 4 GeV 2 , the initial quark distribution remains mostly unchanged 
whether there are zero gluons or, conversely, the gluons are unconstrained. 

The cross sections resulting from fit A (in column 5 of Table ^j) are about 25% of the 
evolved BI cross sections. The diffractive fractions obtained from this fit, using the CTEQ3M 
entries in Table [I], are 2.9% (2.5%) for W ( Z ) production, as shown in Table [|. 

In fit D, there is an enormous amount of glue initially, and this significantly affects the 
evolution of the quarks in the pomeron from = 4 GeV 2 to the vector boson mass. The 
result is an increase in the cross section, as can be seen in column 6 of Table || This happens 
even though the initial quark distribution is smaller than in fit A. Even so, the cross sections 
are still smaller, by a factor of about 2, than the ones from evolved BI pomeron parton 
distributions. The rates from fit D are 4.7% (3.8%) for W ( Z ) production. These rates are 
somewhat lower than those obtained by Kunszt and Stirling ]nj, who mostly used quark 
distributions that at large f3 fall less steeply than ours. 

The data [f|| from which our fits were extracted used a conservative cut on the pomeron 
momentum, x™ ax = 0.01. The pomeron flux factor allows for the xp dependence, but to 
ensure maximum compatibility with the ZEUS data without the assumption of standard 
Regge behavior, the same cut should be applied to the cross sections in hadron-hadron 
collisions. This results in the cross sections in Table [3j, which therefore represent our most 
accurate prediction of diffractive W and Z production, given only the assumption of hard 
scattering factorization, which of course we wish to test. Notice that with this cut the 
diffractive cross sections are over an order of magnitude smaller than with Xp ax = 0.1. The 
percentages obtained with this cut on xp for W ( Z ) production are 0.07% (0.03%) and 0.08% 
(0.03%) for fits A and D, respectively, as shown in Table f|The large reduction is presumably 
due to the fact that we are not far from an effective kinematic limit: the cut on xp gives a 
maximum proton-pomeron energy of 180 GeV, and partons typically do not carry the whole 
of the energy of their parent hadrons. 
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4.2 Why are the fractions smaller than from BI? 

Although the data used in our fits support a “hard” quark distribution^ in the pomeron, we 
predict that the diffractive W and Z cross sections are much smaller, by a factor of 3 to 5, 
than those predicted by Bruni and Ingelman, who also used hard quark distributions. This 
factor arises as an accumulation of several modest factors that all change the ratios in the 
same direction. 

• A factor 0.8 because of the larger inclusive cross sections when one uses CTEQ3M 
instead of the obsolete EHLQ1 distributions in the proton. 

• A factor 0.7 for the effect of the evolution of the parton densities in the pomeron. 

• A factor 0.7 for the use of the Donnachie-Landshoff flux factor instead of the Ingelman- 
Schlein flux factor, when the momentum sum is kept fixed. 

• A factor 0.5 because the data indicate that the quarks give a contribution to the 
momentum sum of 0.5 (with the DL normalization), instead of unity as assumed by 
Bruni and Ingelman. 


4.3 Lepton distributions for W production at the Tevatron 

In this section, we present our results for W production, but now with cuts on the emitted 
lepton l. Specifically, we calculate the electron’s (or positron’s) rapidity (y) distribution 
from Eq. (||) for the diffractive process, and Eq. (|l^) for the inclusive one. For the parton 
distributions in the pomeron, we use our five fits, Eqs. (|l8|) -(^2|), evolved up to the W mass. 

We impose cuts on the lepton that are appropriate for measurements m at the CDF 
experiment. A cut of 20 GeV was imposed on the Et of the emitted lepton so that the 
integration region for Et is 20 GeV < Et < \M W . The parameter xp was integrated over 
the range M^/s < xp < 0.01. 

Fig. |3] shows our results for W~ production. The diffractive cross sections using fits B 
and D differ only by a few percent and are represented by solid curves which overlap, while 
those using fits A and C are denoted by the dashed and dot-dashed curves, respectively. 
These curves exhibit a strong fall-off in the region y e > —0.25 that is a consequence of the 
requirement of a rapidity gap. The overlap between the two solid curves suggests that, at 
least for this particular process in this kinematic region, the extra soft term in fit D for the 
quark in Eq. (|2l|) does not contribute much. The lower dotted curve is the inclusive cross 
section rescaled by 10 -3 . 

The diffractive cross section is about 3% to 4.4% of the nondiffractive one at the leftmost- 
edge of the plots (at y = —3) depending on the fit used. The cross sections using fit B or fit 
D are about 45% larger than those using fits A and C at y = —3, with the ratio decreasing 
as y increases. Since we only include the qq —> W channel in the calculations, and the quark 
densities are nearly the same for the five fits because of the small contribution of the soft 
term, this difference is mostly due to the effect of evolution. The presence of a large amount 
of glue in fits B and D increases the evolved quark densities relative to those in fits A and 


10 This agrees with the HI results 
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Figure 3: Rapidity Distribution of e in W production. A cut Xp < 0.01 was used. 


C. This tendency is even more pronounced for the SG fit (upper dotted curve), where the 
glue is concentrated at large (3. 

The corresponding cross sections for W + are shown in Fig. [|. The cross sections are larger 
than for the W~, because a valence up quark from the proton can be used to make a W + , 
especially at large negative rapidities. In the plot, the rapidity gap exists for y e + > —1.6. 
As in the case of W~ production, the diffractive cross section using fit B or D is larger than 
the one with fit A or C, by about 35% at y e + = — 4, decreasing with y e +. 


4.4 Comparison to CDF data for W production 


The CDF collaboration has presented preliminary data on diffractive W production from pp 
collisions at a/s = 1800 GeV [ill], |12|] . The W's are produced with a rapidity gap in the region 
2.4 < 1771 < 4.2. They find that the fraction of diffractive to non-diffractive W production is 
H R w = [1.15 ± O.bl(stat) ± 0.23 (syst)}%. 

In Table ^]we present our diffractive fractions using Eq. (]5|) with x^ ax = 0.017, which we 
determined from the CDF plot 0 relating r] max and xp. The fractions are about a factor 
of two smaller than the measured rates but are all within the experimental uncertainties. 
They are computed with the diffracted hadron being allowed to be either the proton or the 
antiproton. 
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Fit SG 



Figure 4: Rapidity Distribution of e + in W + production. A cut Xp < 0.01 was used 



Fit A 

Fit B 

Fit C 

Fit D 

W+ + W~ 

0.56% 

0.67% 

0.50% 

0 .66% 


Table 5: Diffractive fractions using CDF cuts. 



















5 Diffractive Jets 


In this section, we present our results for jet production. We imposed the following cuts on 
the jet cross sections. These represent the effect of appropriate experimental cuts 
and of cuts to improve the significance of the signal. 


m 


• We require that two jets are produced in the same half of the detector, i.e., y\U 2 > 0, 
where y t is the rapidity of jet i. This eliminates the region where the jets are in 
opposite hemispheres, since that region is well populated by non-diffractive events but 
is relatively unpopulated by diffractive events, because of the rapidity gap requirement. 

• Each jet is required to have a transverse energy Et greater than 20 GeV. This ensures 
that we are definitely in the perturbative region for the jets, but the cut could be 
relaxed. 

• Each jet’s rapidity satisfies \y\ > y cut = 1.8. 

Next, we integrated over the rapidity of one of the jets to obtain a single jet distribution, 
but still subject to the above cuts on the other jet. Eqs. (||) and O were used for the 
diffractive and inclusive cross sections, respectively, with the parton distributions evolved 
to the scale E % The integration limits used for the Et integral were E^ m = 20 GeV and 
F7p ax = y r s/(e~ y + e y ), while the xp integral was performed up to xp = 0.01. In the following 
discussion, we will denote the rapidity of the final state jet by y^ instead of y. 

The resulting cross sections are shown in Fig. |||. There are no points in the middle part 
of the plot because of the rapidity cut used, as described above. The lower (upper) solid 
curve, which results from using fit D (fit B), is about 10 times larger than both the dashed 
and dot-dashed curves, in which fits A and C, respectively, were used. This reflects the 
sensitivity of this particular type of cross section to the gluon content of the pomeron. The 
lower dotted curve, which is symmetric about y — 0, represents the inclusive cross section 
scaled down by a factor of 10~ 3 . 

We also show in Fig. |S] the cross section obtained when the pomeron parton density with 
a singular gluon is used (fit SG), the upper dotted curve. The resulting cross section is about 
1.5 to 2.5 times larger than that obtained using fit B or D. 

The diffractive jet fractions are shown in Fig. || r where R x 100 is plotted as a function of 
2/jet, with R = ^j e t,inci/^ ]e * • As in Fig. |], the solid curves correspond to the rates when fits B 
and D are used, while the dashed and dot-dashed curves denote the rates for fits A and C, 
respectively. For the non-singular distributions, one finds that the rates R are largest when 
fit B is used, varying from 4.5% to 2.2%. With fit D, which differs from fit B by the absence 
of a soft term in the quark distribution, the rates are about tenths of a percent lower. The 
rates obtained with fits A and C are almost identical and range from 0.43% to about 0.18%. 
The rates are largest at y^ et = —4 then decrease as y^ et increases. Of course, the large rates 
for the distributions B and D, both with the large glue, directly result from the fact that 
there is a gluon induced subprocess. For comparison, we also show the even higher rates 
obtained using the singular gluon fit SG (dotted curve), which vary from 6.7% to 5.2%. 

We end this section by making comparisons with preliminary data on diffractive dijet 
production from CDF and DO at a/s = 1800 GeV. CDF has measured dijet data both with a 
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Figure 5: Rapidity Distribution of Jet Cross Sections. 
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Fit D 



Fit A 
Fit C 
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Figure 6: Diffractive Jet Production Rates. 








Fit A 

Fit B 

Fit C 

Fit D 

g.jet.ditt 

^jetjncl 

0.47% 

5.6% 

0.47% 

5.4% 


Table 6: Diffractive fractions for dijet production using CDF cuts (rapidity gap). 



Fit A 

Fit B 

Fit C 

Fit D 

O-jet.ditt 

^jetjincl 

0.33% 

4.45% 

0.37% 

4.24% 


Table 7: Diffractive fractions for dijet production using CDF cuts (Roman pots). 


rapidity gap requirement [10] and with Roman pots [11] along the antiproton beam direction. 
In the first case, the cross section for dijets produced opposite a rapidity ( 77 ) gap in the region 
2.4 < \rj\ < 4.2 is measured. Each jet is required to have a minimum Et of 20 GeV and 
rapidity \rj\ > 1.8. They also measure the dijet cross section without a rapidity gap, i.e., 
what we refer to in this paper as the inclusive cross section. The diffractive fraction they 
measure is jT2H Rjj = [0.62 ± 0.04 ± 0.09]%. The fractions that we obtain using the above 
cuts are shown in Table [|. The rates obtained with fit D or B are about nine times larger 
while those obtained with fit C or A are about 25% smaller than the measured value. Our 
calculation assumes that either the antiproton or the proton is diffracted. 

With their Roman pot triggered diffractive sample, CDF has measured a diffractive 
fraction of Rjj = [0.109 ± 0.003 ± 0.016]%. The data in this sample correspond to £p in 
the range 0.05 < x-p < 0.1, with the jets having minimum Et of 10 GeV. The fractions we 
obtain using the same kinematic cuts are presented in Table [?j. The ones obtained with fits 
D, B are an order of magnitude larger than the data, while those obtained with fits C, A 
are about 3-4 times larger. Our calculation assumes that only the antiproton is diffracted. 
Since our diffractive parton densities are fitted to data in which the proton may be excited, 
our predictions should be reduced by about (15 ± 10)% H before being compared with data 
in which the isolated proton is detected. However, even after this reduction, our predictions 
are still well above the measured diffractive fraction. 

Finally, DO also has some preliminary data |l4j on diffractive dijet production. They 
require a rapidity gap opposite the dijets which have E™ in = 12 GeV and \rjj et \ > 1.6. The 
diffractive fraction they measure is Rjj = [0.67 ± 0.05]%. Our calculated fractions are 
shown in Table ^]. Again, the fractions obtained using fits D, B are an order of magnitude 
larger than the data, while those obtained with fits C, A are consistent with the data. Our 
calculation assumes that either the antiproton or the proton is diffracted. 



Fit A 

Fit B 

Fit C 

Fit D 

0-jet,ditt 

^.jetqincl 

0.8% 

10.4% 

0.9% 

10.0% 


Table 8: Diffractive fractions for dijet production using DO cuts (rapidity gap). 
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6 Conclusions 


We have presented fits to the DIS diffractive structure functions measured by the ZEUS col¬ 
laboration, together with a lowest order calculation of the resulting diffractive cross sections 
for vector boson production and jet production from pp interactions at the Tevatron. Since 
we used pomeron parton distributions fitted to data on diffractive DIS, the rates represent 
a realistic prediction of the cross sections, given the assumption of factorization. 

The quark distributions are fairly well determined, at least in overall size, and we see that 
diffractive W and Z production is predicted to be a few percent of the inclusive cross section, 
at least if suitable cuts are made to select kinematic configurations that are preferentially 
populated in diffractive events. 

We have derived fractional rates R = cr dlff /cr incl which are much smaller than those 
obtained in the benchmark studies |16| of vector boson production. The rates we have 
calculated are 0.08% (0.03%) for W ( Z) production using fit D, and a cut xp < 0.01. For 
the case of lepton decays of vector bosons, we are able to increase these fractions to several 
percent by requiring the lepton to have large rapidity, on the side opposite the rapidity 
gap. Substantial deviations from these rates would indicate a breakdown of hard scattering 
factorization for diffractive processes. In fact, the preliminary data from the Tevatron appear 
to be consistent with our calculations, within large errors. 

The wide range of gluon distributions that is permitted by the DIS data yields predictions 
for the diffractive jet cross section at the Tevatron (see Fig. ^|) which differ by an order of 
magnitude. Given the assumption of factorization, the smaller jet cross sections, from fits 
A and C, represent the lowest reasonable cross sections since they correspond to a pomeron 
without glue at Q q = 4GeW. However, the fits with the large gluon initial distributions 
(fits B, D, and SG) are preferred for the ZEUS DIS data, in agreement with the results 0 
from the HI collaboration. 

For jet production, rates of 4.5% to 2.1% in the rapidity range —4 < y- ]et < —1.8 were 
calculated using fit D or B (with the large glue). The corresponding rates for the calculations 
using fit A or C (with the small gluon distribution), in the same rapidity range, vary from 
0.43% to 0.18%. We are working on using data 0 on diffractive photoproduction of jets to 
reduce greatly the uncertainty in the gluon density. This will improve the accuracy of our 
predictions, particularly for diffractive jet production. We can already see from the results 
reported by the ZEUS collaboration || that parton distributions with a substantial amount 
of glue are preferred. 

As we saw at the end of Sect. [|, only the lowest of the jet cross sections we compute for 
the Tevatron is consistent with the data. These are the cross sections for which there is no 
initial glue in the parton densities. In the case of the Roman pot data, our calculated cross 
sections are all well above the experimental value; these are the data for which the cuts can 
be most directly implemented in our calculations. These results combined with the ZEUS 
measurement of a substantial gluon component in the pomeron already are suggestive of a 
breakdown of factorization. 

With a modest increase in accuracy, comparison with actual data for diffractive W, Z 
and jet cross sections will permit a good test of the factorization hypothesis. As explained 
in Ref. [[lj , the results of such a test will be important as a probe of the space-time structure 
of hadron-hadron scattering at high energy. 
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